Population exposure to multiple air pollutants and its compound episodes in Europe.

Air pollution remains as a substantial health problem, particularly regarding the combined health risks arising from simultaneous exposure to multiple air pollutants. However, understanding these combined exposure events over long periods has been hindered by sparse and temporally inconsistent monitoring data. Here we analyze daily ambient PM2.5, PM10, NO2 and O3 concentrations at a 0.1-degree resolution during 2003-2019 across 1426 contiguous regions in 35 European countries, representing 543 million people. We find that PM10 levels decline by 2.72% annually, followed by NO2 (2.45%) and PM2.5 (1.72%). In contrast, O3 increase by 0.58% in southern Europe, leading to a surge in unclean air days. Despite air quality advances, 86.3% of Europeans experience at least one compound event day per year, especially for PM2.5-NO2 and PM2.5-O3. We highlight the improvements in air quality control but emphasize the need for targeted measures addressing specific pollutants and their compound events, particularly amidst rising temperatures.

Ambient temperature and risk of cardiovascular and respiratory adverse health outcomes: a nationwide cross-sectional study from Spain.

AIMS: We assessed the association of temperature and temperature variability with cause-specific emergency hospitalizations and mortality from cardiovascular and respiratory diseases in Spain, as well as the effect modification of this association by individual and contextual factors. METHODS AND RESULTS: We collected data on health (hospital admissions and mortality), weather (temperature and relative humidity), and relevant contextual indicators for 48 Spanish provinces during 2004-2019. The statistical analysis was separately performed for the summer (June-September) and winter (December-March) seasons. We first applied a generalized linear regression model with quasi-Poisson distribution to estimate daily province-specific temperature-health associations, and then we fitted multilevel multivariate meta-regression models to the evaluate effect modification of the contextual characteristics on heat- and cold-related risks. High temperature increased the risk of mortality across all cardiovascular and respiratory diseases, with the strongest effect for hypertension (relative risk (RR) at 99th temperature percentile vs. optimum temperature: 1.510 [95% empirical confidence interval {eCI} 1.251 to 1.821]), heart failure (1.528 [1.353 to 1.725]), and pneumonia (2.224 [1.685 to 2.936]). Heat also had an impact on all respiratory hospitalization causes (except asthma), with similar risks between pneumonia (1.288 [1.240 to 1.339]), acute bronchitis and bronchiolitis (1.307 [1.219 to 1.402]), and chronic obstructive pulmonary disease (1.260 [1.158 to 1.372]). We generally found significant risks related to low temperature for all cardiovascular and respiratory causes, with heart failure (RR at 1st temperature percentile vs. optimum temperature: 1.537 [1.329 to 1.779]) and chronic obstructive pulmonary disease (1.885 [1.646 to 2.159]) exhibiting the greatest risk for hospitalization, and acute myocardial infarction (1.860 [1.546 to 2.238]) and pneumonia (1.734 [1.219 to 2.468]) for mortality. Women and the elderly were more vulnerable to heat, while people with secondary education were less susceptible to cold compared to those not achieving this educational stage. Results from meta-regression showed that increasing heating access to the highest current provincial value (i.e. 95.6%) could reduce deaths due to cold by 59.5% (57.2 to 63.5). CONCLUSION: Exposure to low and high temperatures was associated with a greater risk of morbidity and mortality from multiple cardiovascular and respiratory conditions, and heating was the most effective societal adaptive measure to reduce cold-related mortality.

Exposure to low and high temperatures increases the risk of morbidity and mortality from several cardiovascular and respiratory diseases, especially among the elderly. Increasing access to heating could substantially reduce cold-related mortality burden.

Remeasuring the influence of ageing on heat-related mortality in Spain, 1980 to 2018.

Climate change and population ageing are converging challenges that are expected to significantly worsen the health impacts of high temperatures. We aimed to remeasure the implications of ageing for heat-related mortality by comparing time trends based on chronological age (number of years already lived) with those derived from the application of state-of-the-art demographic methodology which better captures the dynamics of evolving longevity: prospective age (number of years still to be lived). We conducted a nationwide time-series analysis of 13 regions in Spain over 1980-2018 using all-cause mortality microdata for people aged 65+ and annual life tables from the Spanish National Institute of Statistics, and daily mean temperatures from E-OBS. Based on confounder-adjusted quasi-Poisson regression with distributed lag non-linear models and multivariate meta-analysis in moving 15-year timeslices, we assessed sex-specific changes in absolute risk and impacts for heat-related mortality at extreme and moderate temperatures, for chronological and prospective age groups. In the conventional chronological age analysis, absolute risk fell over the study period (e.g. females, extreme heat: -54%; moderate heat: -23%); after accounting for rising longevity, the prospective age analysis, however, found a smaller decline in risk for extreme heat (-15%) and a rise for moderate heat (+46%). Additionally, while the chronological age analysis suggested a shift in mortality towards higher ages, the prospective age analysis showed that over the study period, people of largely the same (prospective) age were impacted. Further, the prospective age analysis revealed excess risk in females (compared to males) rose from 20% to 27% for extreme heat, and from 40% to 70% for moderate heat. Assessing the implications of ageing using a prospective age perspective showed the urgency of re-doubling risk reduction efforts, including accelerating healthy ageing programs that incorporate climate considerations. The age patterns of impacts suggested that such actions have the potential to mitigate ageing-related heat-health threats to generate climate change-ready, healthy societies.

The effect of temporal data aggregation to assess the impact of changing temperatures in Europe: an epidemiological modelling study.

Background: Daily time-series regression models are commonly used to estimate the lagged nonlinear relation between temperature and mortality. A major impediment to this type of analysis is the restricted access to daily health records. The use of weekly and monthly data represents a possible solution unexplored to date. Methods: We temporally aggregated daily temperatures and mortality records from 147 contiguous regions in 16 European countries, representing their entire population of over 400 million people. We estimated temperature-lag-mortality relationships by using standard time-series quasi-Poisson regression models applied to daily data, and compared the results with those obtained with different degrees of temporal aggregation. Findings: We observed progressively larger differences in the epidemiological estimates with the degree of temporal data aggregation. The daily data model estimated an annual cold and heat-related mortality of 290,104 (213,745-359,636) and 39,434 (30,782-47,084) deaths, respectively, and the weekly model underestimated these numbers by 8.56% and 21.56%. Importantly, differences were systematically smaller during extreme cold and heat periods, such as the summer of 2003, with an underestimation of only 4.62% in the weekly data model. We applied this framework to infer that the heat-related mortality burden during the year 2022 in Europe may have exceeded the 70,000 deaths. Interpretation: The present work represents a first reference study validating the use of weekly time series as an approximation to the short-term effects of cold and heat on human mortality. This approach can be adopted to complement access-restricted data networks, and facilitate data access for research, translation and policy-making. Funding: The study was supported by the ERC Consolidator Grant EARLY-ADAPT (https://www.early-adapt.eu/), and the ERC Proof-of-Concept Grants HHS-EWS and FORECAST-AIR.

Countrywide analysis of heat- and cold-related mortality trends in the Czech Republic: growing inequalities under recent climate warming.

BACKGROUND: Only little is known about trends in temperature-mortality associations among the most vulnerable subgroups, especially in the areas of central and eastern Europe, which are considered major climatic hotspots in terms of heatwave exposure. Thus, we aimed to assess trends in temperature-related mortality in the Czech Republic by sex, age and cause of death, and to quantify the temporal evolution of possible inequalities. METHODS: We collected daily time series of all-cause (1987-2019) and cause-specific (1994-2019) mortality by sex and age category, and population-weighted daily mean 2-metre temperatures for each region of the Czech Republic. We applied a quasi-Poisson regression model to estimate the trends in region-specific temperature-mortality associations, with distributed lag non-linear models and multivariate random-effects meta-analysis to derive average associations across the country. We then calculated mortality attributable to non-optimal temperatures and implemented the indicator of sex- and age-dependent inequalities. RESULTS: We observed a similar risk of mortality due to cold temperatures for men and women. Conversely, for warm temperatures, a higher risk was observed for women. Results by age showed a clear pattern of increasing risk due to non-optimum temperatures with increasing age category. The relative risk (RR) related to cold was considerably attenuated in most of the studied subgroups during the study period, whereas an increase in the RR associated with heat was seen in the overall population, in women, in the age category 90+ years and with respect to respiratory causes. Moreover, underlying sex- and age-dependent inequalities experienced substantial growth. CONCLUSIONS: Our findings suggest ongoing adaptation to cold temperatures. Mal/adaptation to hot temperatures occurred unequally among population subgroups and resulted in growing inequalities between the sexes and among age categories.

Ambient temperature and seasonal variation in inpatient mortality from respiratory diseases: a retrospective observational study.

Background: The seasonal fluctuation in mortality and hospital admissions from respiratory diseases, with a winter peak and a summer trough, is widely recognized in extratropical countries. However, little is known about the seasonality of inpatient mortality and the role of ambient temperature remains uncertain. We aimed to analyse the association between ambient temperature and in-hospital mortality from respiratory diseases in the provinces of Madrid and Barcelona, Spain. Methods: We used data on daily hospitalisations, weather (ie, temperature and relative humidity) and air pollutants (ie, PM2.5, PM10, NO2 and O3) for the Spanish provinces of Madrid and Barcelona during 2006-2019. We applied a daily time-series quasi-Poisson regression in combination with distributed lag non-linear models (DLNM) to assess, on the one hand, the seasonal variation in fatal hospitalisations and the contribution of ambient temperature, and on the other hand, the day-to-day association between temperature and fatal hospital admissions. The analyses were stratified by sex, age and primary diagnostic of hospitalisation. Findings: The study analysed 1 710 012 emergency hospital admissions for respiratory diseases (mean [SD] age, 60.4 [31.0] years; 44.2% women), from which 103 845 resulted in in-hospital death (81.4 [12.3] years; 45.1%). We found a strong seasonal fluctuation in in-hospital mortality from respiratory diseases. While hospital admissions were higher during the cold season, the maximum incidence of inpatient mortality was during the summer and was strongly related to high temperatures. When analysing the day-to-day association between temperature and in-hospital mortality, we only found an effect for high temperatures. The relative risk (RR) of fatal hospitalisation at the 99th percentile of the distribution of daily temperatures vs the minimum mortality temperature (MMT) was 1.395 (95% eCI: 1.211-1.606) in Madrid and 1.612 (1.379-1.885) in Barcelona. In terms of attributable burden, summer temperatures (June-September) were responsible for 16.2% (8.8-23.3) and 22.3% (15.4-29.2) of overall fatal hospitalisations from respiratory diseases in Madrid and Barcelona, respectively. Women were more vulnerable to heat than men, whereas the results by diagnostic of admission showed heat effects for acute bronchitis and bronchiolitis, pneumonia and respiratory failure. Interpretation: Unless effective adaptation measures are taken in hospital facilities, climate warming could exacerbate the burden of inpatient mortality from respiratory diseases during the warm season. Funding: European Research Council Consolidator Grant EARLY-ADAPT, European Research Council Proof-of-Concept Grants HHS-EWS and FORECAST-AIR.

Drivers of the time-varying heat-cold-mortality association in Spain: A longitudinal observational study.

BACKGROUND: A number of studies have reported reductions in mortality risk due to heat and cold over time. However, questions remain about the drivers of these adaptation processes to ambient temperatures. We aimed to analyse the demographic and socioeconomic drivers of the downward trends in vulnerability to heat- and cold-related mortality observed in Spain during recent decades (1980-2018). METHODS: We collected data on all-cause mortality, temperature and relevant contextual indicators for 48 provinces in mainland Spain and the Balearic Islands between Jan 1, 1980, and Dec 31, 2018. Fourteen contextual indicators were analysed representing ageing, isolation, urbanicity, heating, air conditioning (AC), house antiquity and ownership, education, life expectancy, macroeconomics, socioeconomics, and health investment. The statistical analysis was separately performed for the range of months mostly causing heat- (June-September) and cold- (October-May) related mortality. We first applied a quasi-Poisson generalised linear regression in combination with distributed lag non-linear models (DLNM) to estimate province-specific temperature-mortality associations for different periods, and then we fitted univariable and multivariable multilevel spatiotemporal meta-regression models to evaluate the effect modification of the contextual characteristics on heat- and cold-related mortality risks over time. FINDINGS: The average annual mean temperature has risen at an average rate of 0·36 °C per decade in Spain over 1980-2012, although the increase in temperature has been more pronounced in summer (0·40 °C per decade in June-September) than during the rest of the year (0·33 °C per decade). This warming has been observed, however, in parallel with a progressive reduction in the mortality risk associated to both hot and cold temperatures. We found independent associations for AC with heat-related mortality, and heating with cold-related mortality. AC was responsible for about 28·6% (31·5%) of the decrease in deaths due to heat (extreme heat) between 1989 and 1993 and 2009-2013, and heating for about 38·3% (50·8%) of the reductions in deaths due to cold (extreme cold) temperatures. Ageing (ie, proportion of population over 64 years) attenuated the decrease in cold-related mortality. INTERPRETATION: AC and heating are effective societal adaptive measures to heat and cold temperatures. This evidence holds important implications for climate change health adaptation policies, and for the projections of climate change impacts on human health.

The Direct and Indirect Influences of Interrelated Regional-Level Sociodemographic Factors on Heat-Attributable Mortality in Europe: Insights for Adaptation Strategies.

BACKGROUND: Heat is a significant cause of mortality, but impact patterns are heterogenous. Previous studies assessing such heterogeneity focused exclusively on risk rather than heat-attributable mortality burdens and assume predictors are independent. OBJECTIVES: We assessed how four interrelated regional-level sociodemographic predictors-education, life expectancy, the ratio of older to younger people (aging index), and relative income-influence heterogeneity in heat-attributable mortality burdens in Europe and then derived insights into adaptation strategies. METHODS: We extracted four outcomes from a temperature-mortality study covering 16 European countries: the rate of increase in mortality risk at moderate and extreme temperatures (moderate and extreme slope, respectively), the minimum mortality temperature percentile (MMTP), and the underlying mortality rate. We used structural equation modeling with country-level random effects to quantify the direct and indirect influences of the predictors on the outcomes. RESULTS: Higher levels of education were directly associated with lower heat-related mortality at moderate and extreme temperatures via lower slopes and higher MMTPs. A one standard deviation increase in education was associated with a -0.46±0.14, -0.41±0.12, and 0.41±0.12 standard deviation (±standard error) change in the moderate slope, extreme slope, and MMTP, respectively. However, education had mixed indirect influences via associations with life expectancy, the aging index, and relative income. Higher life expectancy had mixed relations with heat-related mortality, being associated with higher risk at moderate temperatures (0.33±0.11 for the moderate slope; -0.19±0.097 for the MMTP) but lower underlying mortality rates (-0.72±0.097). A higher aging index was associated with higher burdens through higher risk at extreme temperatures (0.13±0.072 for the extreme slope) and higher underlying mortality rates (0.93±0.055). Relative income had relatively small, mixed influences. DISCUSSION: Our novel approach provided insights into actions for reducing the health impacts of heat. First, the results show the interrelations between possible vulnerability-generating mechanisms and suggest future research directions. Second, the findings point to the need for a dual approach to adaptation, with actions that explicitly target heat exposure reduction and actions focused explicitly on the root causes of vulnerability. For the latter, the climate crisis may be leveraged to accelerate ongoing general public health programs. https://doi.org/10.1289/EHP11766.

Heat-related mortality in Europe during the summer of 2022.

Over 70,000 excess deaths occurred in Europe during the summer of 2003. The resulting societal awareness led to the design and implementation of adaptation strategies to protect at-risk populations. We aimed to quantify heat-related mortality burden during the summer of 2022, the hottest season on record in Europe. We analyzed the Eurostat mortality database, which includes 45,184,044 counts of death from 823 contiguous regions in 35 European countries, representing the whole population of over 543 million people. We estimated 61,672 (95% confidence interval (CI) = 37,643-86,807) heat-related deaths in Europe between 30 May and 4 September 2022. Italy (18,010 deaths; 95% CI = 13,793-22,225), Spain (11,324; 95% CI = 7,908-14,880) and Germany (8,173; 95% CI = 5,374-11,018) had the highest summer heat-related mortality numbers, while Italy (295 deaths per million, 95% CI = 226-364), Greece (280, 95% CI = 201-355), Spain (237, 95% CI = 166-312) and Portugal (211, 95% CI = 162-255) had the highest heat-related mortality rates. Relative to population, we estimated 56% more heat-related deaths in women than men, with higher rates in men aged 0-64 (+41%) and 65-79 (+14%) years, and in women aged 80+ years (+27%). Our results call for a reevaluation and strengthening of existing heat surveillance platforms, prevention plans and long-term adaptation strategies.

Daylight saving time affects European mortality patterns.

Daylight saving time (DST) consists in a one-hour advancement of legal time in spring offset by a backward transition of the same magnitude in fall. It creates a minimal circadian misalignment that could disrupt sleep and homoeostasis in susceptible individuals and lead to an increased incidence of pathologies and accidents during the weeks immediately following both transitions. How this shift affects mortality dynamics on a large population scale remains, however, unknown. This study examines the impact of DST on all-cause mortality in 16 European countries for the period 1998-2012. It shows that mortality decreases in spring and increases in fall during the first two weeks following each DST transition. Moreover, the alignment of time data around DST transition dates revealed a septadian mortality pattern (lowest on Sundays, highest on Mondays) that persists all-year round, irrespective of seasonal variations, in men and women aged above 40.

Evolution of temperature-attributable mortality trends looking at social inequalities: An observational case study of urban maladaptation to cold and heat.

BACKGROUND: To date, little is known about the temporal variation of the temperature-mortality association among different demographic and socio-economic groups. The aim of this work is to investigate trends in cold- and heat- attributable mortality risk and burden by sex, age, education, marital status, and number of household occupants in the city of Turin, Italy. METHODS: We collected daily time-series of temperature and mortality counts by demographic and socio-economic groups for the period 1982-2018 in Turin. We applied standard quasi-Poisson regression models to data subsets of 25-year moving subperiods, and we estimated the temperature-mortality associations with distributed lag non-linear models (DLNM). We provided cross-linkages between the evolution of minimum mortality temperatures, relative risks of mortality and temperature-attributable deaths under cold and hot conditions. RESULTS: Our findings highlighted an overall increase in risk trends under cold and heat conditions. All-cause mortality at the 1st percentile increased from 1.15 (95% CI: 1.04; 1.28) in 1982-2006 to 1.24 (95% CI: 1.11; 1.38) in 1994-2018, while at the 99th percentile the risk shifted from 1.51 (95% CI: 1.41; 1.61) to 1.59 (95% CI: 1.49; 1.71). In relation to social differences, women were characterized by greater values in respect to men, and similar estimates were observed among the elderly in respect to the youngest subgroup. Risk trends by educational subgroups were mixed, according to the reference temperature condition. Finally, individuals living in conditions of isolation were characterized by higher risks, with an increasing vulnerability throughout time. CONCLUSIONS: The overall increase in cold- and heat- related mortality risk suggests a maladaptation to ambient temperatures in Turin. Despite alert systems in place increase public awareness and improve the efficiency of existing health services at the local level, they do not necessarily prevent risks in a homogeneous way. Targeted public health responses to cold and heat in Turin are urgently needed to adapt to extreme temperatures due to climate change.

Mortality risk attributable to high and low ambient temperature in Pune city, India: A time series analysis from 2004 to 2012.

BACKGROUND: Exposure to high and low ambient temperatures is associated with morbidity and mortality across the globe. Most of these studies assessing the effects of non-optimum temperatures on health and have been conducted in the developed world, whereas in India, the limited evidence on ambient temperature and health risks and has focused mostly on the effects of heat waves. Here we quantify short term association between all temperatures and mortality in urban Pune, India. METHODS: We applied a time series regression model to derive temperature-mortality associations based on daily mean temperature and all-cause mortality records of Pune city from year January 2004 to December 2012. We estimated high and low temperature-mortality relationships by using standard time series quasi-Poisson regression in conjunction with a distributed lag non-linear model (DLNM). We calculated temperature attributable mortality fractions for total heat and total cold. FINDINGS: The analysis provides estimates of the total mortality burden attributable to ambient temperature. Overall, 6∙5% [95%CI 1.76-11∙43] of deaths registered in the observational period were attributed to non-optimal temperatures, cold effect was greater 5.72% [95%CI 0∙70-10∙06] than heat 0∙84% [0∙35-1∙34]. The gender stratified analysis revealed that the highest burden among men both for heat and cold. CONCLUSION: Non-optimal temperatures are associated with a substantial mortality burden. Our findings could benefit national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately due to climate change.

Projections of temperature-attributable mortality in Europe: a time series analysis of 147 contiguous regions in 16 countries.

BACKGROUND: Europe has emerged as a major climate change hotspot, both in terms of an increase in seasonal averages and climate extremes. Projections of temperature-attributable mortality, however, have not been comprehensively reported for an extensive part of the continent. Therefore, we aim to estimate the future effect of climate change on temperature-attributable mortality across Europe. METHODS: We did a time series analysis study. We derived temperature-mortality associations by collecting daily temperature and all-cause mortality records of both urban and rural areas for the observational period between 1998 and 2012 from 147 regions in 16 European countries. We estimated the location-specific temperature-mortality relationships by using standard time series quasi-Poisson regression in conjunction with a distributed lag non-linear model. These associations were used to transform the daily temperature simulations from the climate models in the historical period (1971-2005) and scenario period (2006-2099) into projections of temperature-attributable mortality. We combined the resulting risk functions with daily time series of future temperatures simulated by four climate models (ie, GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, and MIROC5) under three greenhouse gas emission scenarios (ie, Representative Concentration Pathway [RCP]2.6, RCP6.0, and RCP8.5), providing projections of future mortality attributable fraction due to moderate and extreme cold and heat temperatures. FINDINGS: Overall, 7·17% (95% CI 5·81-8·50) of deaths registered in the observational period were attributed to non-optimal temperatures, cold being more harmful than heat by a factor of ten (6·51% [95% CI 5·14-7·80] vs 0·65% [0·40-0·89]), and with large regional differences across countries-eg, ranging from 4·85% (95% CI 3·75-6·00) in Germany to 9·87% (8·53-11·19) in Italy. The projection of temperature anomalies by RCP scenario depicts a progressive increase in temperatures, more exacerbated in the high-emission scenario RCP8.5 (4·54°C by 2070-2099) than in RCP6.0 (2·89°C) and RCP2.6 (1·67°C). This increase in temperatures was transformed into attributable fraction. Projections consistently indicated that the increase in heat attributable fraction will start to exceed the reduction of cold attributable fraction in the second half of the 21st century, especially in the Mediterranean and in the higher emission scenarios. The comparison between scenarios highlighted the important role of mitigation, given that the total attributable fraction will only remain stable in RCP2.6, whereas the total attributable fraction will rapidly start to increase in RCP6.0 by the end of the century and in RCP8.5 already by the middle of the century. INTERPRETATION: The increase in heat attributable fraction will start to exceed the reduction of cold attributable fraction in the second half of the 21st century. This finding highlights the importance of implementing mitigation policies. These measures would be especially beneficial in the Mediterranean, where the high vulnerability to heat will lead to an imbalance between the decreasing cold and increasing heat-attributable mortality. FUNDING: None.

Trade-offs between short-term mortality attributable to NO2 and O3 changes during the COVID-19 lockdown across major Spanish cities.

The emergence of the COVID-19 pandemic forced most countries to put in place lockdown measures to slow down the transmission of the virus. These lockdowns have led to temporal improvements in air quality. Here, we evaluate the changes in NO2 and O3 levels along with the associated impact upon premature mortality during the COVID-19 lockdown and deconfinement periods along the first epidemic wave across the provincial capital cities of Spain. We first quantify the change in pollutants solely due to the lockdown as the difference between business-as-usual (BAU) pollution levels, estimated with a machine learning-based meteorological normalization technique, and observed concentrations. Second, instead of using exposure-response functions between the pollutants and mortality reported in the literature, we fit conditional quasi-Poisson regression models to estimate city-specific associations between daily pollutant levels and non-accidental mortality during the period 2010-2018. Significant relative risk values are observed at lag 1 for NO2 (1.0047 [95% CI: 1.0014 to 1.0081]) and at lag 0 for O3 (1.0039 [1.0013 to 1.0065]). On average NO2 changed by -51% (intercity range -65.7 to -30.9%) and -36.4% (-53.7 to -11.6%), and O3 by -1.1% (-20.2 to 23.8%) and 0.6% (-12.4 to 23.0%), during the lockdown (57 days) and deconfinement (42 days) periods, respectively. We obtain a reduction in attributable mortality associated with NO2 changes of -119 (95% CI: -273 to -24) deaths over the lockdown, and of -53 (-114 to -10) deaths over the deconfinement. This was partially compensated by an increase in the attributable number of deaths, 14 (-72 to 99) during the lockdown, and 8 (-27 to 50) during the deconfinement, associated with the rise in O3 levels in the most populous cities during the analysed period, despite the overall small average reductions. Our study shows that the potential trade-offs between multiple air pollutants should be taken into account when evaluating the health impacts of environmental exposures.

Correction: Heat-related mortality trends under recent climate warming in Spain: A 36-year observational study.

[This corrects the article DOI: 10.1371/journal.pmed.1002617.].

Reversal of the seasonality of temperature-attributable mortality from respiratory diseases in Spain.

A growing number of epidemiological studies have recently assessed temporal variations in vulnerability and/or mortality attributable to hot and cold temperatures. However, the eventual changes in the seasonal distribution of temperature-attributable mortality remain unexplored. Here, we analyse countrywide daily time-series of temperature and mortality counts from respiratory diseases by sex, age group and province of residence during the period 1980-2016 in Spain. We show the complete reversal of the seasonality of temperature-attributable mortality, with a significant shift of the maximum monthly incidence from winter to summer, and the minimum monthly incidence from early and late summer to winter. The reversal in the seasonal distribution of the attributable deaths is not driven by the observed warming in both winter and summer temperatures, but rather by the very large decrease in the risk of death due to cold temperatures and the relatively much smaller reduction due to hot temperatures. We conclude that the projected decrease in the number of moderate and extreme cold days due to climate warming will not contribute to a further reduction of cold-attributable respiratory deaths.